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ABSTRACT 


The  depth  of  the  marine  atmospheric  boundary  layer  (MABL)  is  an  important 
parameter  for  both  scientific  and  operational  meteorological  applications.  The  depth  of 
the  marine  boundary  layer  has  a  significant  influence  on  the  atmospheric  dynamics  in  the 
coastal  zone.  Knowledge  of  the  depth  of  stratocumulus-topped  boundary  layers  (STBLs) 
will  enable  coastal  operations  to  more  accurately  anticipate  weather,  and  electromagnetic 
propagation  conditions.  This  study  develops  a  satellite  remote  sensing  technique  for 
determining  the  height  of  MABLs  topped  with  stratocumulus  clouds. 

Validation  of  the  technique  using  coastal  rawinsonde  dataset  from  the  Monterey 
Area  Ship  Track  (MAST)  experiment  revealed  that  an  assumption  of  41%  cloud  with  a 
moist  lapse  rate  equal  to  -7.0°C/km  had  the  best  overall  fit  to  the  data.  However,  for 
shallow  boundary  layers  with  depths  below  400m  the  most  accurate  assumption  was  75% 
cloud  with  a  moist  lapse  rate  equal  to  -6.5°C/km.  The  application  of  this  technique  to 
sounding  data  returned  an  overall  BL  depth  accuracy  of  50m  while  the  satellite 
application  returned  an  overall  accuracy  of  65m.  A  sensitivity  analysis  of  both  surface 
and  cloud-top  temperature  revealed  that  a  1/2°C  change  in  either  temperature  resulted  in 
an  error  of  60-70m  in  boundary  layer  depth. 
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I.  INTRODUCTION 

The  depth  of  the  marine  atmospheric  boundary  layer  (MABL)  is  an  important 
parameter  for  both  scientific  and  operational  meteorological  applications.  The  depth  of 
the  marine  boundary  layer  has  a  significant  influence  on' the  atmospheric  dynamics  in  the 
coastal  zone.  Since  the  MABL  is  often  topped  by  stratocumulus  clouds,  knowledge  of 
the  depth  of  stratocumulus-topped  boundary  layers  (STBLs)  will  enable  coastal 
operations  to  more  accurately  anticipate  weather,  and  electromagnetic  propagation 
conditions.  A  satellite  remote  sensing  technique  for  determining  the  height  of  MABLs 
topped  with  stratocumulus  clouds  is  developed  in  this  thesis. 

A.  APPLICATION  OF  BOUNDARY  LAYER  DEPTH  DATA 

The  value  of  knowing  boundary  layer  height  can  be  seen  when  considering 
propagation  conditions  of  surface  based  radar  platforms.  Even  without  complete 
knowledge  of  ducting  profile  shapes  in  the  region  knowledge  of  boundary  layer  height  is 
important  using  the  assumption  that  the  boundary  layer  is  well  mixed.  Boundary  layer 
(BL)  height  is  directly  related  to  the  base  of  the  trapping  layer  associated  with 
atmospheric  ducts  that  can  enhance  the  propagation  of  electromagnetic  (EM)  energy.  The 
base  of  the  trapping  layer  is  referred  to  as  the  “optimum-coupling  layer”.  When 
conditions  are  such  that  the  atmosphere  acts  as  a  wave-guide  or  duct,  the  propagation  of 
EM  waves  is  confined  to  a  narrow  region.  Operational  sensors  that  detect  radar  waves 
and  other  EM  energy  can  operate  more  effectively  when  favorable  ducting  conditions  are 
present. 
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Surface  vessels  can  experience  enhanced  EM  detection  conditions  in  the  presence 
of  surface-based  ducts.  The  “optimum-coupling  layer”,  which  is  equivalent  to  the  height 
of  the  well-mixed  portion  of  the  MABL,  is  related  to  the  type  of  duct  present.  Figure  1-1 
illustrates  the  height  difference  for  a  surface-based  duct  (Zj)  and  for  an  elevated  duct  (Z2). 
The  surface-based  duct  normally  has  a  lower  optimum-coupling  layer  (Davidson,  2000). 
The  contrast  between  the  two  types  of  ducts  and  their  impact  on  EM  propagation 
conditions  can  be  seen  in  Figure  1-2  which  shows  greater  detection  ranges  for  surface- 
based  ducts  vice  elevated  ducts  (Rogers,  L.T.,  1999). 


Different  Types  of  Ducts 


Figure  1-1:  Examples  of  elevated  and  surface-based  ducting 
conditions. 
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Different  types  of  ducting 

Evaporation  ducting  Surface-based  ducting 
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Figure  1-2:  Radar  image  of  elevated  and  surface-based  ducting  conditions. 

While  boundary  layer  depth  information  is  of  critical  importance  for  determining 
the  type  and  effect  of  ducting  conditions,  the  moisture  content  above  the  BL  must 
also  be  know  in  order  to  determine  if  a  trapping  layer  occurs,  i.e.,  to  determine  the 
shape  of  the  refractivity  (M)  profile.  Due  to  the  presence  of  strong  subsidence  the 
air  above  the  BL  is  generally  very  dry.  Figure  1-3  shows  a  sounding  from  the 
Arabian  Gulf,  which  clearly  shows  the  dry  nature  of  the  atmosphere  above  the 
boundary  layer  in  the  presence  of  subsidence  (Davidson,  2000).  Ducting  conditions 
are  very  common  especially  over  eastern  ocean  basins  where  strong  subsidence 
from  a  semi-permanent  high  is  common.  Figure  1-4  shows  the  importance  of  the 
knowledge  and  application  of  boundary  layer  depth  to  the  determination  of  ducting 
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conditions  can  be  seen  when  considering  the  prevalence  of  ducting  conditions 
around  the  world.  (Ortenburger  et  al,  1985) 


Airborne/Atmosphere  Characterization  Test 
T  and  Td  Rocketsonde  profile  :  1  sec  rate 


Mean  fall  rate  =  1.7  m/sec 


Figure  1-3  :  Example  of  the  nature  of  moisture  profiles 
in  and  above  the  BL  ducting  in  a  region  of  strong 
subsidence. 


Rosenthal  et  al,  [1997]  showed  that  it  is  possible  to  obtain  a  statistical  correlation 
of  duct  height  (BL  height)  with  cloud-top  temperature  (TCT)  off  the  coast  of 
California  to  within  approximately  133m  of  accuracy  80%  of  the  time  (Rosenthal, 
J.S.,  2000).  The  correlation  of  Tct  to  duct  height  implies  that  a  remote  sensing 


4 


technique  for  determining  BL  depth  under  strong  subsidence  inversions  would 
provide  accurate,  real-time  estimation  of  ducting  conditions. 


Figure  1-4:  World  wide  frequency  of  ducting  occurrence,  (from 
Ortenburger,  L.N.,  et  al,  1985) 

B.  METHODS  OF  CLOUD  HEIGHT  DETERMINATION 

Direct  measurements  of  boundary  layer  height  are  often  difficult  and  impractical  to 
make  and  are  therefore  relatively  rare.  This  has  led  to  the  development  of  several  methods 
by  which  the  height  of  the  MABL  may  be  directly  and  indirectly  determined  in  cloudy  and 
in  cloud  free  atmospheres.  In  the  past  techniques  for  determining  the  height  of  the  MABL 
were  limited  to  land  and  ship  launched  radiosondes.  However,  several  techniques  have 
been  developed  to  accomplish  the  same  goal  in  a  more  efficient  manner.  Some  examples 
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of  these  techniques  include  ground-based  SOund  Detection  And  Ranging  (SODAR) 
(Gaynor  and  Mandics,  1978)  and  Light  Detection  and  Ranging  (LIDAR)  (Eloranta  et  al., 
1975).  These  techniques  require  multiple  stations  and  even  in  optimal  circumstances  they 
fail  to  provide  a  complete  view  of  the  marine  boundary  layer  due  to  their  limited  spatial 
extent. 

The  development  of  satellite-derived  retrieval  methods  provides  an  expanded 
observation  area  and  therefore  a  more  complete  picture  of  boundary  layer  depth.  Many 
satellite  techniques  exist  for  measuring  boundary  layer  properties  including  multispectral 
techniques  for  cloud  free  boundary  layers  (Kren,  1987).  This  study  will  focus  on  cloudy 
boundary  layers  since,  as  Kuze  and  Chance  [1994]  have  shown,  clouds  cover  more  than 
50%  of  the  Earth’s  surface. 

There  are  several  types  of  cloud  height  retrieval  methods  that  utilize  satellite  data. 
These  techniques  can  be  broken  up  into  two  general  categories.  The  first  group  extracts 
cloud  top  height  in  a  direct  way  from  satellite  measurements  including  stereoscopic 
analysis  and  cloud  shadow  measurements.  Stereoscopic  analysis  uses  offset  image  pairs 
that  are  correlated  to  give  a  horizontal  shift  that  is  turned  into  a  cloud  top  height. 
However,  this  technique  is  only  applicable  when  near-coincident  satellite  image  pairs  are 
available  for  comparison  and  is  therefore  restricted  to  geostationary  imagery  (Shenk  et  al., 
1975).  Cloud  shadow  measurements  use  the  given  solar  and  satellite  geometry  to  calculate 
cloud  height  from  the  length  and  direction  of  observed  cloud  shadows.  This  technique 
applies  only  to  clouds  with  limited  vertical  (Hasler  et  al.,  1991). 
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The  second  category  uses  the  physical  characteristics  of  the  atmosphere  to  model 
cloud  height.  Some  of  these  techniques  include  the  use  of  atmospheric  absorption  bands 
(C02  or  02),  and  bispectral  estimations.  Atmospheric  absorption  techniques  use  model 
optical  path  length  and  therefore  height  from  within  an  absorbing  region  of  the 
electromagnetic  spectrum.  Typical  meteorological  satellites  such  as  the  Advanced  Very 
High  Resolution  Radiometer  (AVHRR)  instrument  on  the  NOAA  polar  orbiters  employs 
neither  the  correct  wavelength  regions  nor  the  narrow  bandwidths  needed  to  estimate 
cloud  top  height.  Bispectral  methods  on  the  other  hand  derive  cloud-top  height  by 
comparing  cloud-top  brightness  temperature  (TB)  with  satellite  derived  vertical  soundings 
such  as  from  the  Tiros  Operational  Vertical  Sounder  (TOVS).  Unfortunately,  satellite 
soundings  lack  the  vertical  resolution  that  is  necessary  for  accurate  estimations  of  STBL 
cloud  height.  Another  option  is  to  compare  TB  with  sounding  output  from  numerical 
weather  model  temperature  fields.  These  types  of  fields  nominally  have  temperature  errors 
of  1-2°C  and  this  type  of  application  generally  yields  cloud  height  errors  of  300-500 
meters  (Simpson,  et  al,  2000). 

The  use  of  satellite  data  for  broad  spatial  measurements  of  the  depth  of  the  STBL 
has  the  potential  to  significantly  improve  the  understanding  of  large-scale  STBL  properties 
and  dynamics.  This  study  will  present  and  test  a  technique  for  estimating  the  depth  of  the 
STBL  by  using  a  vertical  thermodynamic  model  that  calculates  boundary  layer  depth  from 
independent  measurements  of  sea  surface  temperature  (Ts)  and  satellite  derived  cloud-top 
brightness  temperature. 
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In  Chapter  II  of  this  study  the  theoretical  background  behind  the  structure  and 
dynamics  of  the  STBL  as  well  as  the  theory  behind  satellite  derived  cloud-top  brightness 
temperatures  will  be  presented.  Chapter  III  will  present  the  datasets  utilized  in  this  study 
as  well  as  the  specifics  of  the  retrieval  technique.  Chapter  FV  will  present  the  pertinent 
assumptions  needed  to  accurately  estimate  the  depth  of  the  STBL  and  the  results  of  the 
technique  applied  to  atmospheric  soundings  and  satellite  derived  TB.  Chapter  V  covers  the 
conclusions  drawn  from  this  thesis  research  and  recommendations  for  further  study. 


II.  BACKGROUND 


With  the  proper  assumptions,  it  is  possible  to  indirectly  determine  the  depth  of  the 
STBL  using  sea  surface  temperature  and  satellite  derived  cloud-top  brightness 
temperature.  The  technique  requires  brightness  temperatures  from  AVHRR  channel  4 
(11  pm  wavelength),  and  an  assumption  of  vertical  cloud  fraction  in  order  to  determine 
the  height  of  transition  from  the  dry  to  moist  (pseudo-adiabatic)  lapse  rate  within  the 
STBL.  This  chapter  covers  the  structure  and  dynamics  of  cloudy  marine  boundary  layers, 
and  the  basic  radiative  transfer  process  that  determines  the  measurement  of  cloud-top 
brightness  temperature. 

A.  LAPSE  RATES  IN  THE  STBL 

The  Earth’s  atmosphere  cools  with  increasing  height  due  to  the  fact  that  it  is 
heated  primarily  from  below  by  longwave  radiation  emitted  from  the  Earth.  The  rate  of 
cooling  is  known  as  the  lapse  rate  of  the  atmosphere.  The  lapse  rate  of  the  atmosphere  is 
associated  with  the  mathematical  slope  of  the  line  of  temperature  vs.  altitude.  By 
convention  a  decrease  of  temperature  for  increasing  altitude  is  denoted  as  a  positive  lapse 
rate. 

The  technique  described  here  assumes  adiabatic  lapse  rates  within  the  STBL. 
However,  the  atmosphere  near  the  Earth’s  surface  is  adiabatic  only  under  certain 
conditions.  These  conditions  include  heavy  cloud  cover,  which  limits  radiative  heating  of 
the  surface,  and  turbulence  sufficient  to  mix  the  air  and  smooth  out  temperature  gradients. 
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In  other  words,  the  atmosphere  may  not  be  adiabatic  due  to  insufficient  mixing  and/or 
solar  heating  of  the  surface. 

These  two  processes  produce  highly  variable  temperature  behavior  in  the  lowest 
layers  of  the  atmosphere.  Solar  heating  of  the  surface  is  most  pronounced  over  land- 
covered  areas.  Over  the  ocean  the  high  heat  capacity  of  water  prevents  significant  diurnal 
warming  and  cooling  of  the  ocean  surface  even  during  cloud  free  days.  Therefore,  the 
STBL  can  be  said  to  meet  these  two  criteria  most  of  the  time,  and  an  adiabatic  lapse  rate 
is  therefore  a  good  assumption. 

1.  Dry  Adiabatic  Lapse  Rate 

The  dry  adiabatic  lapse  rate  is  the  rate  at  which  air  temperature  decreases  with 
increasing  altitude  for  an  air  parcel  that  is  unsaturated.  The  decrease  in  temperature  is 
due  to  the  adiabatic  expansion  (occurring  without  loss  or  gain  of  heat)  of  a  parcel  of  air. 
An  expanding  parcel  of  air  will  cool  as  the  volume  increases  due  to  conservation  of 
energy.  The  value  for  the  dry  adiabatic  lapse  rate  is  approximately  -9.8°C/km. 

2.  Moist  Adiabatic  Lapse  Rate 

The  moist  adiabatic  lapse  rate  is  the  loss  of  temperature  with  increasing  altitude 
for  an  air  parcel  that  is  saturated.  If  a  saturated  air  parcel  is  adiabatically  expanded,  the 
temperature  will  decrease,  but  this  decrease  causes  the  parcel  to  become  increasingly 
supersaturated.  Saturation  leads  to  condensation,  and  the  latent  heat  released  in  this 
process  warms  the  parcel.  The  warming  from  the  latent  heat  release  partially  offsets  the 
loss  of  temperature  due  to  adiabatic  expansion  leading  to  a  lapse  rate  that  is  smaller  in 
magnitude  than  the  dry  adiabatic  lapse  rate. 
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It  can  be  seen  in  Figure  2-1  that  the  pseudo-adiabatic  lapse  rate  is  a  nonlinear 
function  due  to  the  addition  of  latent  heat.  However,  this  lapse  rate  is  essentially  constant 
below  1.5km.  Since  the  STBL  rarely  exceeds  1.5km  in  height,  a  constant  value  is  used 
for  the  pseudo-adiabatic  lapse  rate  in  this  procedure.  (Hsieh,  1987) 


Figure  2-1 :  Atmospheric  lapse  rates, 
(from  Hsieh,  1987) 


B.  STRUCTURE  OF  THE  STBL 

For  cloudless  boundary  layers  the  majority  of  the  turbulence  is  produced  at  the 
surface  by  thermal  convection  and/or  vertical  wind  shear.  Surface  turbulence  leads  to  the 
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upward  transport  of  momentum  and  moisture  in  the  layer.  In  the  STBL  an  additional 
source  of  turbulence  is  added  by  the  presence  of  the  cloud.  Stratiform  clouds  at  the  top  of 
the  marine  boundary  layer  result  in  dynamics  that  are  significantly  more  complicated  than 
for  the  cloud  free  boundary  layer.  The  STBL  is  produced  and  maintained  due  to  the 
balance  of  several  processes.  The  nature  of  these  processes  will  be  discussed  in  the 
following  sections. 

1.  Cloud-Top  Inversion 

STBLs  form  under  inversion  layers  that  are  due  to  the  descending,  or  subsiding  air 
in  the  atmosphere  associated  with  high-pressure  areas.  The  inversion  layer  itself  is  a 
product  of  the  adiabatic  heating  of  the  descending  air  and  the  layer  just  above  the 
subsidence  inversion  is  usually  dry.  Within  an  inversion  layer  there  is  a  reversal  of  the 
normal  lapse  rate  where  temperature  increases  with  altitude.  The  STBL  is  most  prevalent 
beneath  strong  subsidence  inversions  over  subtropical  oceans.  These  areas  include  eastern 
subtropical  ocean  basins  and  mid  to  high  latitude  oceans  where  marine  stratiform  clouds 
are  common.  Specifically,  the  area  in  the  North  Pacific  Ocean  off  the  California  coast  is 
typically  covered  with  marine  stratocumulus  clouds  during  the  summer  months  when  the 
region  is  dominated  by  a  persistent  sub-tropical  high  pressure  system. 

The  degree  of  turbulent  mixing  within  the  boundary  layer  determines  the 
distribution  of  moisture,  which  influences  the  thickness  of  the  cloud  and  can  change  the 
overall  depth  of  the  boundary  layer.  The  following  is  a  breakdown  of  the  processes  that 
produce  the  mixing  required  to  maintain  stratocumulus  clouds.  (Stull,  1988) 
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Surface-based  Free  Convection:  Free  convection  occurs  when  cold  air  is 
advected  over  a  warmer  surface  since  the  condition  of  cold  air  over  warm  water 
results  in  an  unstable  atmosphere.  The  maintenance  of  stratocumulus  clouds  and 
fog  off  the  coast  of  southern  California  is  strongly  influenced  by  surface-based 
free  convection. 

Shear-generated  Mechanical  Turbulence:  Shear  turbulence  is  associated  with 
strong  BL  winds  that  can  cause  convection  (mixing)  that  is  sufficiently  strong 
enough  to  mix  the  moisture  from  the  ocean  surface  to  the  cloud  layer.  The 
efficient  redistribution  of  surface  properties  caused  by  mixing  is  an  important 
factor  in  the  maintenance  of  marine  stratocumulus  clouds. 

Cloud-top  Radiative  Cooling:  The  radiative  heat  loss  at  the  top  of 
stratocumulus  clouds  is  a  significant  factor  in  the  maintenance  of  the  STBL.  The 
cooling  at  the  tops  of  stratocumulus  clouds  creates  pockets  of  cold  air  that  sink 
and  mix  with  lower  levels  and  help  maintain  the  condensation  necessary  to 
produce  cloud  droplets.  Cloud-top  cooling  is  a  process  that  is  important  to  the 
maintenance  of  stratocumulus  clouds  as  long  as  there  are  not  higher  clouds  that 
can  reduce  the  cooling  of  the  lower  cloud  deck  by  radiating  energy  down  onto  the 
cloud  top. 

Cloud-base  Radiative  Heating:  The  process  of  radiative  heating,  by  radiation 
emitted  from  the  surface  and  absorbed  at  the  cloud-base,  plays  an  important  role 
in  the  destabilization  of  the  cloud  layer  while  stabilizing  the  subcloud  layer. 
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However,  it  is  generally  much  weaker  than  cloud-top  cooling  due  to  the  smaller 
temperature  difference  between  layers. 

Cloud-layer  Shear:  Shear  generated  in  the  cloud-layer  generally  produces 
mechanical  turbulence  and  mixing  near  the  top  of  the  cloud.  The  turbulence 
resulting  from  cloud-layer  shear  produces  small  scale  entrainment  mixing  that  can 
decouple  the  cloud  unless  there  are  other  turbulent  processes  present  to  more 
evenly  distribute  moisture  and  temperature  throughout  the  BL. 

2.  Cloud  Base  or  LCL 

The  altitude  at  which  saturation  occurs  inside  an  air  parcel,  as  it  rises,  expands  and 
cools,  is  called  the  lifting  condensation  level  (LCL).  The  LCL  corresponds  to  the  base  of 
the  clouds  in  the  STBL.  The  value  of  the  LCL  is  not  constant  in  the  STBL  and  its 
variation  is  governed  by  the  degree  of  mixing  or  turbulence  within  the  boundary  layer  and 
by  a  destructive  process  known  as  decoupling. 

Decoupling  is  the  process  by  which  stratocumulus  clouds  become  disconnected 
from  the  lower  regions  of  the  marine  boundary  layer.  When  the  boundary  layer  is 
decoupled  a  small  inversion  develops  below  the  cloud  base,  which  inhibits  mixing  from 
the  ocean  surface  to  the  cloud  and  creates  a  lower,  moist  layer  beneath  a  drier,  cloud 
mixed  layer.  Decoupling  occurs  as  a  result  of  two  primary  processes  in  the  boundary  layer 
(Tjemstrom,  1998). 

a.  Decoupling  by  Solar  Absorption  (Diurnal  Decoupling) 

The  process  of  diurnal  decoupling  is  caused  by  variations  in  mixing  due  to 
cloud-top  radiative  cooling.  During  both  day  and  night  the  longwave  radiation  that  is 
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emitted  upward  from  the  cloud  top  is  generally  higher  than  the  incoming  energy  received 
from  the  dry  above  the  cloud.  The  net  loss  of  radiative  energy  causes  an  imbalance  that 
cools  the  top  of  the  cloud  layer  and  leads  to  thermally  induced  instability  since  the  top  of 
the  cloud  is  now  cooler  than  the  bottom.  The  instability  is  enhanced  under  daytime 
heating  as  shortwave  solar  radiation  penetrates  into  the  cloud  which  warms  it  from 
within. 

The  effect  of  solar  radiation  can  cause  a  secondary  inversion  to  form  at  the 
cloud  base  if  the  surface  flux  is  weak.  The  secondary  inversion  separates  the  cloud  from 
the  surface  and  a  second  turbulent  layer  forms  in  the  STBL.  The  diurnal  fluctuation  of 
cloud  thickness  due  to  this  process  plays  an  important  role  in  the  understanding  of  the 
dynamics  of  the  STBL. 

b.  Decoupling  Due  to  Precipitation 

The  second  decoupling  process  is  due  to  precipitation  out  the  bottom  of 
the  cloud.  As  described  above,  vertical  cooling  rate  of  the  atmosphere  is  modified  by  the 
release  of  latent  heat  as  water  vapor  condenses.  However,  as  drizzle  falls  out  of  the  cloud 
some  of  it  evaporates  in  the  air  below  leading  to  a  cooler,  moister  sub-cloud  layer.  The 
development  of  a  sub-cloud  layer  causes  a  stable  layer  to  form  near  the  cloud  base  that 
prevents  the  boundary  layer  from  being  well  mixed. 

3.  The  Surface  (Sea  and  Air  Temperature) 

The  surface  air  temperature  (TA)  and  sea  surface  temperature  (Ts)  are  usually 
close  in  magnitude.  If  a  gradient  in  temperature  occurs,  fluxes  develop  which  produce 
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differences  in  TA  and  Ts  such  that  they  seldom  precisely  equal  in  value.  There  are  several 
physical  process  which  can  affect  the  differences  in  TA  and  Ts.  These  processes  include: 
Convection:  Convection  is  a  process  by  which  heat  is  transferred  due  to  the 
physical  movement  of  air.  Convective  processes  help  to  minimize  temperature 
fluctuations  by  vertically  mixing  the  air  near  the  surface. 

Conduction:  The  transfer  of  heat  through  matter  by  communication  of  kinetic 
energy  from  particle  to  particle  is  known  as  conduction.  Air  is  made  up  of  many 
dispersed  molecules  that  inhibit  the  process  of  conduction  making  air  an  effective 
insulator.  The  process  of  conduction  is  negligible  when  compared  to  the  other 
process  in  the  atmosphere,  which  are  much  larger  in  magnitude. 

Radiation:  Radiation  is  the  transfer  of  heat  via  electro-magnetic  energy  which  is 
transferred  to  the  atmosphere  through  absorption  by  gaseous  constituents  such  as 
water  vapor  and  carbon  dioxide.  Radiation  from  the  ocean  has  a  very  strong 
moderating  effect  on  the  air  temperature  over  the  ocean’s  surface,  which  is  due  to 
the  large  heat  capacity  of  the  ocean.  The  high  heat  capacity  of  the  ocean  enables  it 
to  gain  and  lose  energy  without  large  changes  in  Ts- 

Advection:  In  many  cases,  advection  is  the  most  influential  of  the  four  processes 
in  amplifying  the  sea/air  temperature  difference.  When  an  air  mass  moves  out  of 
its  generation  area  and  moves  over  ocean  waters  with  different  properties,  strong 
fluxes  of  temperature  and  moisture  result.  These  fluxes  act  to  reduce  the 
difference  between  the  sea  surface  and  the  surface  air  temperature. 
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Klein  and  Hartmann  [1993]  showed  the  annual,  climatological  cycle  of  Ts  and  TA 
difference  to  be  very  small  in  the  Californian  stratus  region.  Figure  2-2  shows  the 
difference  between  TA  and  Ts  is  typically  between  0.5°C  and  1°C.  This  was  done  when 
they  compared  the  seasonal  variation  of  stratus  cloud  amounts  off  the  coast  of  California 
using  shipboard  data.  It  is  therefore  assumed  that  the  difference  between  the  air  and  sea 
surface  temperature  is  nominally  within  1°C  at  any  given  time.  This  small  (1°C) 
temperature  difference  is  primarily  due  to  climatologically  weak  surface  fluxes  of 
sensible  heat.  For  simplicity  in  this  study,  it  is  initially  assumed  here  that  the  sea  surface 
temperature  is  equal  to  the  surface  air  temperature. 


Figure  2-2:  The  annual  cycle  of  Ts  and  TA  difference. 
(Klein  and  Hartmann,  1993) 
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C.  SATELLITE  DETERMINED  CLOUD-TOP  TEMPERATURE 

The  technique  presented  in  this  thesis  uses  the  AVHRR  sensor  to  determine 

cloud-top  temperature.  AVHRR  channel  4  measures  radiance  at  11pm  wavelength. 
Channel  4  is  ideal  for  measuring  surface  or  cloud-top  radiance  due  to  the  lack  of  gaseous 
absorption  in  the  atmosphere  at  this  wavelength.  (Kidder  and  Vonder  Haar,  1995) 

1.  Radiative  Transfer  to  Determine  Cloud-Top  Temperature 
It  is  assumed  that  the  cloud  layer  being  sensed  is  thick  and  that  the  atmosphere 
above  the  cloud  top  is  dry.  From  Schwartzchild’s  equation  the  radiance  at  the  top  of  the 
atmosphere  is: 


L  (X, 0, cp)  =  es (X, Q)B(X, Ts)xd (X)  +  £  B(^, T(p))  — (>:-P)  dp 


dp 


(2-1) 


where  Lt  is  the  spectral  radiance  at  the  top  of  the  atmosphere  as  a  function  of  wavelength 
(A.)  and  direction  (0,cp).  £s(A.,0)B(X,Ts)xd(A.)  is  the  emitted  radiance  from  the  surface  as  a 
function  of  wavelength  and  surface  temperature.  xd(A.)  is  the  total  direct  transmittance 
throughout  the  full  extent  of  the  atmosphere  as  a  function  of  wavelength  and  xd(/-,p)  is  the 
direct  transmittance  from  pressure  p  to  the  top  of  the  atmosphere  (p=0).  B(X,T(p))  is  the 
Planck  blackbody  emittance  as  a  function  of  wavelength,  and  temperature  at  a  given 
pressure.  The  derivative  of  the  direct  transmittance  (xd)  with  respect  to  pressure  is  a 
weighting  function  that  determines  the  contribution  to  Lt  due  to  at  each  vertical  position 
in  the  atmosphere. 
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For  the  STBL  [dxd(A,p)/dp]dp  =  1  over  a  few  meter  thick  layer  thick  near  the 
cloud  top.  Therefore,  nearly  all  of  the  radiance  emitted  from  the  cloud  comes  from  the 
upper  portion  only.  The  surface  term  [ss(L,Ts)  xd(A,)]  goes  to  zero  since  the  cloud  is 
opaque  at  these  wavelengths  and  it  absorbs  all  of  the  energy  from  the  surface.  The  final 
result  is: 

L=B(X,Tct)  (2-2) 

Here  Lt  is  measured  by  the  sensor  at  X  and  the  value  of  TCt  can  be  readily  determined 
from  inversion  of  the  Planck  blackbody  relation. 

2.  Errors  in  Satellite  Determined  TCt  of  the  STBL 

Several  conditions  can  potentially  interfere  with  the  determination  of  TCT  from 
satellite  Tb. 

a.  High  Clouds 

The  presence  of  higher  clouds  (e.g.,  cirrus)  above  the  STBL  produces 
lower  values  of  TCT  due  to  the  lower  emitted  radiance  of  clouds  that  exist  in  the  higher, 
colder  levels  of  the  atmosphere.  The  presence  of  cirrus  clouds  in  particular  can  be  a 
problem  due  to  the  fact  that  cirrus  clouds  are  often  thin  and  difficult  to  detect  without 
multispectral  techniques.  Therefore,  the  presence  of  cirrus  clouds  in  the  vicinity  of  the 
STBL  can  introduce  significant  errors  in  low  level,  satellite  derived  cloud-top 
temperatures. 
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b.  Limb  Darkening 

When  a  measurement  is  taken  from  a  satellite  at  the  periphery  of  its  field 
of  view  the  radiance  that  is  detected  must  pass  through  a  larger  volume  of  atmosphere 
than  radiance  measured  from  nadir.  The  effect  of  greater  atmospheric  interaction  reduces 
the  amount  of  cloud-top  radiance  that  reaches  the  satellite  sensor.  The  AVHRR  sensor 
has  a  scan  angle  from  nadir  of  +/-55.30  (Kidder  and  Vonder  Haar,  1995).  Measurements 
made  at  an  angle  of  55.3°  are  not  only  made  through  an  increased  atmospheric  volume 
but  radiance  emitted  from  the  cloud-top  has  an  angular  dependence  as  well.  The  radiance 
emitted  from  a  cloud  decreases  as  the  view  angle  increases.  Therefore  measurements  of 
Tct  made  at  the  edge  of  a  satellite  pass  are  often  significantly  reduced  in  radiance. 

D.  THE  VALIDITY  OF  A  WELL-MIXED  STBL 

Within  the  STBL  the  potential  temperature  and  water  vapor  mixing  ratio  are 
homogeneous  due  to  turbulent  mixing  within  the  layer.  The  eddies  which  produce  this 
mixing  are  due  to  the  presence  of  buoyancy  and  shear  effects.  The  magnitude  of 
buoyancy  production  within  the  STBL  is  small  due  to  the  high  heat  capacity  of  the  ocean 
which  prevents  large  diurnal  temperature  fluctuations.  Nighttime  radiative  cooling  of  the 
sea  surface  prevents  the  establishment  of  a  near  surface  stable  layer.  Therefore,  the  STBL 
is  often  neutrally  stable  with  only  slight  diurnal  variations.  (Kren,  1987) 

The  weakly  unstable  or  neutral  STBL  is  influenced  by  wind  shear  as  well.  The 
effect  of  wind  shear  on  the  STBL  is  to  produce  turbulent  eddies  which  mix  the  STBL. 
Turbulence  produced  by  both  buoyancy  and  shear  is  generally  sufficient  to  mix  the  STBL. 
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Some  regions  under  synoptically  induced  subsidence  inversions  have  a  stable 
STBL  where  mixing  is  often  suppressed.  These  areas  have  relatively  cold  sea  surface 
temperatures  that  often  comes  from  the  coastal  upwelling  found  along  the  west  coast  of 
continents.  In  these  areas  where  the  STBL  is  not  well  mixed,  using  a  dry  adiabatic  lapse 
rate  will  add  uncertainty  to  the  final  value  of  boundary  layer  depth.  The  accuracy  of  the 
technique  used  in  this  thesis  under  these  conditions  will  be  described  in  the  next  few 
chapters  of  this  text. 

The  technique  presented  in  this  thesis  has  certain  limitations  which  constrain  its 
applicability.  The  technique  requires  knowledge  of  the  temperature  of  both  the  sea 
surface  and  the  cloud-top.  However,  satellite  techniques  for  obtaining  sea-surface 
temperature  require  a  cloud-free  atmosphere,  and  techniques  for  obtaining  cloud-top 
temperature  requires  either  an  atmosphere  with  little  water  vapor  above  the  STBL,  or 
knowledge  of  the  water  vapor  content  so  its  effects  can  be  corrected.  The  technique 
presented  here  has  potential  application  in  regions  of  strong  subsidence  inversions  where 
cloudy  boundary  layers  persist.  Sea  surface  temperature  data  can  be  provided  by  ship 
observations,  recent  satellite  analyses  of  SST  from  clear  sky  periods  or  predictive  ocean 
models. 
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III.  DATA  AND  METHODOLOGY 

A.  MONTEREY  AREA  SHIP  TRACK  (MAST)  EXPERIMENT  DATASET 

The  primary  data  set  for  this  thesis  comes  from  the  Monterey  Area  Ship  Track 
Experiment  (MAST)  (Durkee  et  al,  2000).  MAST  was  a  multi-agency  effort  designed  to 
determine  the  physical  nature  of  ship  tracks  via  the  collection  of  sounding,  aircraft  and 
satellite  data  off  the  coast  of  Central  California.  Figure  3-1  depicts  the  MAST  operating 
area.  This  thesis  will  focus  on  AVHRR  satellite  data  and  radiosonde  data  collected  by  the 
R/V  Glorita  during  the  MAST  experiment. 


Figure  3-1:  MAST  Experiment  1994  Operations  Area. 
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1. 


Radiosonde  Data 


The  R/V  Glorita  launched  an  average  of  six  radiosondes  a  day  at  intervals  of 
roughly  four  hours  for  a  total  of  94  radiosonde  soundings  (designated  GL01-GL94). 
Parameters  measured  with  the  radiosondes  were  pressure,  temperature  and  relative 
humidity.  These  measured  quantities  where  then  used  to  derive  values  for  height, 
dewpoint  and  mixing  ratio.  The  techniques  used  to  derive  these  quantities  are  included  in 
Appendix  B. 

2.  Surface  Observation  Data 

Personnel  onboard  the  R/V  Glorita  made  hourly  measurements  of  sea  surface 
temperature  (Ts),  air  temperature  (TA),  relative  humidity,  as  well  as  wind  speed  and 
direction.  However,  gaps  in  Ts  measurements  did  occur  when  the  temperature  sensor  was 
not  in  the  water.  Additionally  these  measurements  are  amplified  by  general  comments  on 
sky  conditions  made  during  radiosonde  launch  times. 

3.  Satellite  (AVHRR)  Data 

Imagery  from  AVHRR  Channels  1-5  of  the  NOAA-9, 10,  1 1,  and  12  satellites  was 
recorded  during  the  experiment.  Channel  4  (11pm)  brightness  temperatures  from  this 
imagery  were  used  to  determine  the  cloud  top  temperature  at  the  location  of  a  given 
sounding  from  the  R/V  Glorita. 

B.  ANALYSIS  PROCEDURE 

In  order  to  determine  the  depth  of  the  STBL  a  vertical,  thermodynamic  model  was 
constructed.  The  model  assumes  that  unsaturated  air  follows  the  dry  adiabatic  lapse  rate 
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(Td)  until  the  LCL  is  reached.  From  the  LCL  to  the  cloud  top  the  air  is  saturated  and 
follows  a  pseudo  or  moist  adiabatic  lapse  rate  (Tm).  With  a  given  assumption  of  cloud 
thickness  these  two  lapse  rates  can  be  combined  in  order  to  give  a  more  complete 
thermodynamic  picture  of  the  STBL. 

The  combination  of  Td  and  rm  determines  a  cloud-top  height  (Z)  for  a  given 
surface  temperature  (Ts)  and  cloud-top  temperature(TCT).  The  constant  lapse  rate 
equation  T  =  (TCt-Ts)/Z  is  simply  rearranged  algebraically  to  yield  the  formula  for  height 
Z  =  (Tct-Ts)/T .  Given  a  formula  for  height  based  on  temperature  and  lapse  rate  only,  it  is 
now  possible  to  compute  heights  based  on  the  moist  and  dry  adiabatic  lapse  rates  for 
given  values  of  Ter  and  Ts. 

An  assumption  of  vertical  cloud  fraction  must  be  made  in  order  to  give  a  proper 
physical  parameterization  of  how  much  of  the  boundary  layer  follows  Tm  and  how  much 
follows  Td.  Since  the  actual  height  of  the  boundary  layer  is  unknown  an  initial  value,  or 
first  guess,  for  the  boundary  layer  height  must  be  used  in  order  to  determine  the  LCL  or 
cloud  base  height.  The  LCL  in  this  technique  is  derived  from  an  assumed  fraction  of  the 
total,  dry  adiabatic  height  (Zd). 

Figure  3-2  shows  an  example  where  it  is  assumed  that  the  upper  one-third  of  the 
boundary  layer  is  cloud  and  the  remaining  two-thirds  of  the  boundary  layer  is  cloud-free 
air.  The  LCL  is  determined  by  multiplying  the  total  height  of  the  cloud-free  boundary 
layer  (obtained  via  the  dry  adiabatic  lapse  rate)  by  two-thirds.  The  remainder  of  the 
boundary  layer  height  is  calculated  from  the  LCL  by  using  the  pseudo-adiabatic  lapse  rate 
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until  the  value  of  Tct  is  reached.  However,  by  virtue  of  the  fact  that  the  pseudo-adiabatic 
lapse  rate  is  steeper  than  the  dry  lapse  rate  this  process  produces  a  slightly  larger  fraction 
of  cloud  than  is  initially  assumed.  Therefore,  when  an  assumption  of  l/3rd  cloud  is  made 
based  on  the  height  of  a  totally  dry  BL,  the  increase  in  height  from  the  moist-adiabatic 
lapse  rate  produces  a  final  BL  distribution  of  41%  cloud  and  59%  cloud-free  atmosphere. 


nc* 
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Figure  3-2:  The  application  of  the  assumption  of  a  cloud  thickness  of 
41%  and  how  it  is  derived  from  the  dry  adiabatic  height  Zd. 

In  order  to  calculate  the  height  of  the  lower,  cloud-free  portion  of  the  boundary 
layer  the  height  of  the  entire  boundary  layer  (i.e.,  from  Ts  to  TCt)  must  be  computed 
along  the  dry  adiabat  with  the  relation  Zd  =  (TCt-Ts)/  Td  where  Zd  is  the  height  of  the 
boundary  layer  if  it  were  entirely  cloud  free.  The  height  of  the  cloud  free  portion  of  the 
boundary  layer  is  then  simply  the  height  of  the  cloud-free  BL  multiplied  by  a  slight 
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overestimate  of  the  initial  fraction  of  cloud-free  STBL  or  ZCF  =  Zd*(FCF).  Where  ZCF 
denotes  the  initial  estimate  of  the  height  of  the  cloud-free  STBL  and  Fcf  is  the  initial 
estimate  of  the  cloud-free  fraction  of  the  STBL. 

Once  the  initial  height  of  the  cloud  free  portion  of  the  boundary  layer  (which  is 
actually  the  height  of  the  cloud  base)  is  determined,  the  corresponding  temperature  at  that 
height  must  be  known  in  order  to  calculate  the  height  of  the  saturated  portion  of  the 
boundary  layer.  The  temperature  at  the  cloud  base  or  LCL  is  computed  by  the  relation 
Tcf  =  ZCF*rd  +  Ts  which  then  allows  for  the  same  fundamental  lapse  rate  relation  to  be 
used  in  the  upper  two  thirds  of  the  boundary  layer.  This  final  height  is  combined  with  the 
height  of  the  below  cloud  BL  to  yield  the  total  height  of  the  STBL  via  the  following 
relation:  Zbl  =  ZcF  +  (Tct-TcF)/  rm. 

In  order  to  determine  the  accuracy  of  the  technique,  values  of  TA  and  TCt  were 
taken  directly  from  the  Glorita  soundings  and  used  to  determine  the  height  of  the  STBL. 
To  determine  the  accuracy  of  the  technique  using  satellite  data,  the  inversion  temperature 
from  the  sounding  is  replaced  by  the  cloud-top  brightness  temperature  from  the  AVHRR 
channel  4  imagery.  A  total  of  46  soundings  were  used  from  the  MAST  dataset  and  48 
soundings  were  dismissed  outright  from  the  MAST  dataset  due  to  insufficient  MABL 
clouds.  See  Table  3  in  Appendix  A  for  more  detail  on  which  soundings  were  dismissed 
from  the  dataset. 
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IV.  RESULTS 


A.  ASSUMPTIONS 

Two  primary  assumptions  must  be  made  in  order  to  model  the  general  state  of  the 
STBL.  The  amount  of  cloud  must  be  specified  and  the  lapse  rates  of  the  atmosphere  in 
the  cloud  and  below  cloud  must  be  given.  In  the  STBL  it  is  assumed  that  all  lapse  rates 
are  of  constant  value  and  near  adiabatic  in  nature.  The  below  cloud  atmosphere  is 
assumed  to  be  exactly  dry  adiabatic  with  a  lapse  rate  of  Td  =  -9.8°C/km  and  within  cloud 
the  value  of  the  pesudo-adiabatic  lapse  rate  is  initially  given  the  standard  atmospheric 
value  of  rm  =  -6.5°C/km. 

It  was  initially  assumed  that  the  most  common  ratio  of  saturated  to  unsaturated  air 
in  the  STBL  was  2/3rd  cloud  to  l/3rd  unsaturated  (below  cloud)  air.  However,  results 
have  shown  that  a  more  accurate  depiction  of  the  STBL  requires  an  assumption  of  75% 
cloud  for  shallow  boundary  layers  (below  400m  deep)  and  a  value  of  41%  cloud  for 
boundary  layers  deeper  than  400m.  Above  400m  the  value  for  the  pseudo-adiabatic  lapse 
rate  was  modified  to  a  value  of  Tm  =  -7.0°C/km  to  partially  compensate  for  its  non-linear 
nature  as  altitude  increases.  The  reasons  these  assumptions  were  chosen  will  be 
discussed  below. 

B.  SOUNDING  APPLICATION 

In  order  to  validate  the  physical  reasoning  of  the  technique  described  in  Chapter 
IE,  the  method  was  applied  to  the  soundings  from  the  R/V  Glorita.  The  method  was  run 
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with  several  different  values  for  rm  and  cloud  thickness  in  order  to  determine  the  best  set 
of  parameter  values.  The  results  for  each  set  of  assumptions  were  organized  by  the 
difference  between  the  actual  BL  height  and  the  calculated  height  of  the  STBL  (denoted 
as  AH).  Throughout  the  verification  process  the  assumptions  made  about  the  STBL  were 
modified  in  order  to  minimize  the  value  of  AH 

1.  Actual  Depth  Method 

The  distribution  of  AH  with  respect  to  the  actual  BL  depth  was  plotted  in  Figure 
4-1,  which  shows  how  the  two  sets  of  assumptions  (75%  and  41%  cloud)  begin  to  diverge 
in  the  vicinity  of  400  meters.  Based  on  these  criteria  it  was  determined  from  the  sounding 
results  that  there  was  a  discontinuity  in  the  accuracy  of  a  single  set  of  assumptions  at 
around  the  400m  point.  Boundary  layer  cases  with  depth  less  than  400m  are  more 
accurately  modeled  by  an  assumption  of  75%  cloud  with  a  lapse  rate  of  Tm  =  -6.5°C/km. 
These  values  were  the  only  set  of  assumptions  that  returned  a  value  for  AH  <  100m  in  all 
cases  below  400m.  For  boundary  layers  deeper  than  400m  it  was  determined  that  an 
assumption  of  41%  cloud  combined  with  a  lapse  rate  of  Tm  =  -7.0°C/km  returned  the 
most  accurate  height.  These  assumptions  for  soundings  above  400m  produced  a  height 
differential  (AH)  that  was  consistently  below  100m  for  the  majority  of  the  soundings. 

Figure  4-2  shows  how  the  two  sets  of  assumptions  were  combined  into  a  two- 
layer  model  in  order  to  give  a  more  accurate  picture  of  the  total  STBL.  The  height 
produced  by  this  two-layer  procedure  were  plotted  against  the  actual  BL  height  with  the 
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Figure  4-1 :  Scatter  plot  of  the  two  primary 
sets  of  assumptions  applied  all  soundings. 

The  results  of  this  regression  produced  a  line  with  a  slope  of  1 .0283  where  a  slope 
of  1.0  is  a  perfect  correlation.  The  line  was  offset  from  the  origin  by  a  value  of  -41m 
while  the  overall  error  of  the  estimate  of  this  technique  is  56m.  Table  1  in  Appendix  A 
shows  how  several  of  the  different  assumption  sets  affect  the  determination  of  boundary 
layer  height. 
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Figure  4-2:  Linear  regression  of  actual  vs.  calculated 
boundary  layer  height  using  sounding  data.  The  solid  line 
is  the  regression  line  and  N  is  the  number  of  data  points. 

2.  First  Guess  Method 

While  this  process  produces  favorable  results  it  requires  previous  knowledge  of 
the  actual  depth  of  the  STBL.  Since  the  actual  depth  of  the  STBL  is  the  goal  of  this  thesis 
it  cannot  be  used  experimentally  to  determine  the  dividing  line  between  the  two  sets  of 
assumptions.  Therefore,  a  “first  guess”  must  be  made  in  order  to  determine  which  set  of 
assumptions  to  use. 

It  can  be  seen  in  Table  1  that  for  a  single  set  of  assumptions  the  most  accurate 
parameters  over  the  entire  STBL  are  a  41%  cloud  layer  and  pseudo-adiabatic  lapse  rate  of 
Tm  =  -7.0°C/km.  These  assumptions  were  applied  to  the  entire  STBL  and  once  an  initial 
value  of  BL  depth  was  obtained,  the  “shallow”  set  of  assumptions  were  applied  to  the 
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soundings  below  the  depth  of  400m.  Figure  4-3  shows  the  results  from  this  two-step 
process  where  soundings  were  plotted  against  the  actual  BL  depth  and  a  linear  regression 
was  performed.  The  regression  was  found  to  be  even  more  accurate  than  the  original 
application  of  the  assumptions  with  a  slope  of  1.0080,  an  offset  of -24m  and  an  overall 
error  estimate  of  50m. 


Figure  4-3:  Linear  regression  of  actual  vs.  calculated 
boundary  layer  height  using  sounding  data.  The  solid  line  is 
the  regression  line  and  N  is  the  number  of  data  points. 

*  -  Indicates  the  composite  BL  height  was  determined  by 
the  two-step  process  of  the  application  of  BL  assumptions. 

The  soundings  from  the  R/V  Glorita  indicate  that  shallow  BLs,  (ZBL  <  400m) 
have  a  larger  fraction  of  the  BL  that  is  saturated  (cloudy)  air  than  for  deeper  BLs. 
However,  the  height  value  of  400m  is  by  no  means  an  absolute  dividing  line  between  BLs 
with  differing  percentages  of  cloud.  The  MAST  dataset  contains  several  cases  where 
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soundings  with  BL  depths  above  400m  contain  less  cloud  than  75%.  When  the 
assumption  of  41%  cloud  is  applied  to  these  cases  a  significant  underestimate  is  returned. 
This  occurs  due  to  the  fact  that  when  the  assumption  of  41%  cloud  is  applied  to  these 
cases,  the  predominance  of  the  dry  lapse  rate  produces  a  significant  underestimate  in  BL 
depth.  When  these  soundings  are  reevaluated  with  the  shallow  assumption  set  (75% 
cloud),  a  more  accurate  determination  of  the  depth  of  the  STBL  results  due  to  the  steeper 
pseudo-adiabatic  lapse  rate.  Therefore,  the  process  of  making  a  first  guess  at  the  depth  of 
the  BL  transforms  the  dividing  point  for  the  two  sets  of  assumptions  into  a  broader  range 
of  depths.  For  this  dataset  the  effect  was  to  extend  the  depth  from  400m  up  to  460m  in 
cases  where  the  BL  depth  was  severely  underestimated  by  the  assumption  of  41%  cloud. 

3.  Sounding  Case  Studies 

Table  3  lists  all  R/V  Glorita  soundings  along  with  position,  launch  time  and 
general  comments  on  their  analysis.  A  total  of  46  soundings  were  studied  in  the  sounding 
application  part  of  this  thesis  in  order  to  validate  the  various  assumptions  made  about  the 
properties  of  the  STBL.  The  following  is  a  breakdown  of  several  soundings  that  are 
representative  of  three  primary  depth  categories.  These  cases  are  presented  in  detail  in 
order  to  demonstrate  how  the  assumptions  described  above  apply  to  the  various  boundary 
layer  depths  encountered  during  this  study. 

a.  GL16 

Figure  4-4  shows  a  plot  of  sounding  GL16  that  was  launched  at  1 153  UTC 
on  11  June  94  at  36.15°N-122.57°W.  The  surface  pressure  for  this  sounding  is  1015.0 
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mb,  relative  humidity  is  98.8%,  and  winds  are  from  the  NNW  at  0.3m/s  and  increasing  in 
strength  with  height.  Ts  was  not  recorded  for  this  case  and  surface  air  temperature  is 


13.2°C.  The  depth  of  the  boundary  layer  in  this  case  is  330  meters  with  a  cloud  fraction 
of  approximately  75%. 
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Figure  4-4:  Sounding  GL16  with  lapse  rate  assumptions  of 
75%  cloud  and  Tm  =  -6.5°C/km  applied. 

When  the  procedure  is  run  on  the  sounding  GL16  the  first  guess 
assumption  (41%  cloud)  returns  a  BL  height  of  254m  which  corresponds  to  a  AH  of  - 
76m  or  an  underestimate  of  76m.  This  falls  in  the  <400m  category  and  the  technique 
requires  a  shift  in  the  assumption  set  to  the  75%  cloud  assumption.  This  combination  of 
lapse  rates  gives  a  much  more  favorable  result  where  the  calculated  height  is  now  300m 
with  a  AH  of  only  -30m. 
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b. 


GL67 


Figure  4-5  shows  a  plot  of  Sounding  GL67  that  was  launched  at  1148 
UTC  on  22  June  94  at  35.78°N-123.18°W.  The  surface  pressure  for  this  sounding  is 
1019.8  mb,  relative  humidity  is  77.6%,  and  winds  are  from  the  NW  at  9.6m/s.  Surface  air 
temperature  is  14.1°C,  but  Ts  was  not  recorded  for  this  case.  The  depth  of  the  boundary 
layer  in  this  case  is  504  meters. 

When  the  procedure  is  applied  to  sounding  GL67  the  first  guess 
assumption  of  41%  cloud  returns  a  height  value  that  is  well  above  the  400m  height  and 
therefore  these  values  will  be  retained  in  dealing  with  this  sounding.  The  41%  cloud 
assumption  returned  a  calculated  height  of  507m  with  a  AH  of +3m. 

It  is  important  to  note  that  when  this  sounding  is  compared  to  a  satellite 
image  that  is  described  later  it  can  be  seen  that  it  passes  through  an  area  of  reduced 
cloudiness  that  is  immediately  adjacent  to  a  broad  area  of  very  even  marine  stratocumulus 
clouds.  Due  to  the  uniform  nature  of  the  surrounding  cloudiness  it  is  assumed  that  the 
boundary  layer  properties  associated  with  the  surrounding  clouds  is  represented  in  the 
sounding  GL67  even  though  the  sounding  does  not  pass  directly  through  these  clouds. 
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Figure  4-5:  Sounding  GL67  with  lapse  rate  assumptions  of 
41%  cloud  and  rm  =  -7.0°C/km  applied. 

c.  GL68 

Figure  4-6  shows  a  plot  of  Sounding  GL68  that  was  launched  at  1453  UTC 
on  22  June  94  at  35.70°N-123.46°W.  The  surface  pressure  for  this  sounding  is  1021.1 
mb,  relative  humidity  is  84.4%,  and  winds  are  from  the  NW  at  9.0m/s  and  decreasing  in 
strength  with  height.  The  surface  air/sea  temperature  difference  is  TA-TS  =  -1.1°C,  and 
the  surface  air  temperature  is  14.0°C.  The  depth  of  the  boundary  layer  in  this  case  is 
671m  with  a  cloud  fraction  of  approximately  37%. 
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Figure  4-6:  Sounding  GL68  with  lapse  rate  assumptions 
of  41%  cloud  and  rm  =  -7.0°C/km  applied. 

When  the  lapse  rate  procedure  is  applied  to  this  sounding  the  first  guess 
assumption  is  utilized  for  this  case  due  to  the  returned  boundary  layer  depth  of  665m, 
which  corresponds  to  a  AH  of  -6m.  As  a  test  case  the  actual  Ts  value  (Ts  =  15.1°C)  was 
used  in  this  method  and  the  result  was  an  increase  in  the  calculated  depth  of  the  boundary 
layer.  This  initial  temperature  modified  the  result  such  that  an  assumption  of  41%  cloud 
returns  a  calculated  height  of  789m  with  a  AH  of  +118m.  This  indicates  that  a  1°C 
change  in  surface  temperature  modifies  the  height  differential  value  by  approximately 
120m.  The  effect  of  surface  and  cloud-top  temperature  variations  on  the  entire  sounding 
dataset  will  be  discussed  in  section  C  of  this  chapter. 
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d.  GL17 

Figure  4-7  shows  a  plot  of  Sounding  GL17  that  was  launched  at  1448  UTC 
on  11  June  94  at  36.23°N-122.48°W.  The  surface  pressure  for  this  sounding  is  1015.5 
mb,  relative  humidity  is  99.3%,  and  winds  are  from  the  West  at  l.Om/s. 


File  a  gif  7 .  jxt  Percertf  Error  k  -1 B  016S 


Figure  4-7:  Sounding  GL17  with  lapse  rate  assumptions  of 
75%  cloud  and  Tm  =  -6.5°C/km  applied. 

The  surface  air  temperature  for  this  sounding  is  12.2°C  and  the  value  of  Ts  was  not 
recorded  for  this  sounding.  The  depth  of  the  boundary  layer  in  this  case  is  340m  with  a 
cloud  fraction  of  approximately  88%. 

Sounding  GL17  is  relatively  shallow  and  when  the  first  guess  set  of 
assumptions  are  applied  it  returns  a  height  of  242m  and  therefore  the  75%  cloud 
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assumption  is  applied  to  this  sounding.  The  assumption  of  75%  cloud  gives  a  calculated 
BL  height  of  286m  which  corresponds  to  an  underestimate  of  AH  =  -54m. 

e.  GL29 

Figure  4-8  shows  a  plot  of  Sounding  GL29  that  was  launched  at  2347 
UTC  on  12  June  94  at  37.05°N-123.25°W.  The  surface  pressure  for  this  sounding  is 
1016.8  mb,  relative  humidity  is  91.9%,  and  winds  are  from  the  NW  at  7.0m/s  and 
increasing  in  strength  with  height.  The  surface  air  temperature  for  this  sounding  is  14.0°C 
and  the  air/sea  temperature  difference  is  TA-TS  =  -0.6°C.  The  depth  of  the  boundary  layer 
in  this  case  is  377  meters  with  a  cloud  fraction  of  approximately  23.5%. 

This  sounding  initially  meets  the  conditions  necessary  to  be  reevaluated 
with  the  assumptions  for  BLs  below  400m.  The  first  guess  routine  returns  a  BL  depth  of 
277m  with  a  corresponding  height  differential  of  AH  =  -100m.  When  the  assumption  of 
75%  cloud  is  applied  to  this  sounding  is  returns  a  boundary  layer  height  of  327m  which 
corresponds  to  a  height  differential  of  AH  =  -50m. 
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file  =  gl29  t.xt  Petcem  E*ror  =  -13.2364 


Figure  4-8:  Sounding  GL29  with  lapse  rate  assumptions  of  75%  cloud 
and  fm  =  -6.5°C/km  applied. 


C.  TEMPERATURE  VARIATIONS 

It  was  shown  in  the  previous  section  that  variations  in  surface  temperature  can 
have  a  significant  effect  on  the  accuracy  of  the  technique  presented  in  this  thesis.  Figure 
4-9  shows  an  individual  example  where  sounding  GL29  was  modified  such  that  the  lapse 
rate  is  started  at  the  actual  sea  surface  temperature  (Ts  =  14.6°C). 
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Temp  &  Dew  Point  (C) 


Figure  4-9:  Sounding  GL29  corrected  to  start  at  the  actual  SST  of  14.6°C. 

The  steep,  light  blue  lines  represent  the  two  different  values  for  the 
pseudo-adiabatic  lapse  rate. 

A  sensitivity  test  was  made  for  the  entire  sounding  dataset  where  both  surface  and 
cloud-top  temperature  were  modified.  The  results  of  these  modifications  are  compared  to 
the  41%  cloud  assumption  as  well  as  a  “no  cloud”  and  an  “all  cloud”  state.  Figures  4-10 
thru  4-13  demonstrate  how  an  offset  in  surface  and  cloud- top  temperature  (Ts,  Tct)  affects 
the  accuracy  of  the  model.  Figures  4-10  and  4-12  indicate  that  a  %°C  decrease  in  Ts  or  a 
‘/2°C  increase  in  Tct  produces  a  resulting  accuracy  that  is  equivalent  to  that  of  a  model  that 
relies  solely  on  a  dry  adiabatic  lapse  rate  for  the  entire  BL.  That  is,  these  temperature 
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changes  approximate  that  of  a  model  that  utilizes  a  dry  adiabatic  lapse  rate  over  the  entire 
depth  of  the  BL.  Figures  4-1 1  and  4-12  show  how  a  1°C  increase  in  Ts  or  a  1°C  decrease 
in  Tct  produces  a  result  that  is  equivalent  to  that  of  a  model  that  relies  solely  on  a  dry 
adiabatic  lapse  rate  for  the  entire  BL,  i.e.,  a  boundary  layer  completely  filled  with  cloud. 

The  blue  (41%  cloud)  line  as  with  all  of  the  lines  represents  a  single  set  of 
assumptions  applied  to  the  entire  dataset.  The  slight  increase  in  error  or  tail  that  lies  away 
from  and  below  the  main  peak  represents  the  shallow  cases  where  ZBl  <400m.  These 
figures  give  a  sense  of  the  impact  of  inaccuracies  in  both  Ts  and  Tct  where  even  relatively 
small  errors  in  temperature  can  have  a  significant  impact  on  the  accuracy  of  the  technique. 
These  results  imply  that  if  errors  Tcx  or  Ts  are  greater  than  1°C  then  there  is  essentially  no 
advantage  to  breaking  up  the  lapse  rates  when  determining  BL  depth. 
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Figure  4-10:  The  effect  of  a  decrease  in  surface  temperature  compared 
to  a  cloud-free  boundary  layer. 
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Figure  4-11:  The  effect  of  an  increase  in  surface  temperature 
compared  to  an  all-cloud  boundary  layer. 
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Figure  4-12:  The  effect  of  an  increase  in  cloud-top  temperature 
compared  to  a  cloud-free-boundary  layer. 
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Figure  4-13:  The  effect  of  a  decrease  in  cloud- top  temperature 
compared  to  cloud-free  boundary  layer. 
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D.  SATELLITE  APPLICATION 

The  sounding  profiles  from  the  Glorita  were  paired  with  images  from  NOAA- 
9,10,11  and  12  that  were  closest  in  time  to  the  soundings.  Due  to  the  relatively 
homogeneous  nature  of  environment  that  produces  marine  stratocumulus  clouds  a 
maximum  time  difference  of  three  hours  was  chosen  as  the  criteria  to  identify 
sounding/image  pairs. 

Channel  4  of  each  satellite  image  was  analyzed  in  the  region  that  the  sounding 
was  taken  in  order  to  determine  the  cloud  top  temperature.  However,  the  presence  of 
cirrus  clouds  in  the  image  caused  some  cases  to  be  dismissed.  Cirrus  clouds,  being  much 
colder  than  MABL  clouds,  contaminated  the  brightness  temperature  and  produce 
erroneous  results  that  significantly  overestimate  the  depth  of  the  STBL.  Therefore,  this 
method  is  not  valid  for  images  with  significant  cirrus  contamination.  For  detailed 
information  on  which  soundings  were  utilized  for  this  procedure  and  which  were 
dismissed  see  Table  4  in  Appendix  A. 

The  technique  applied  to  the  sounding  dataset  was  applied  to  satellite  imagery  in 
order  to  evaluate  the  accuracy  of  the  application  of  Tb  measurements  that  are  utilized  as 
cloud-top  temperature.  The  value  of  Ts  in  this  part  of  the  study  was  once  again  obtained 
from  the  sounding  in  order  to  better  understand  the  accuracy  of  substituting  TB  for  TCt. 
The  returned  values  of  BL  height  were  plotted  against  the  actual  BL  height  with  the  actual 
BL  height  as  the  independent  variable. 

Figures  4-14  and  4-15  show  a  linear  regression  was  performed  on  the  data  points 
in  order  to  determine  the  accuracy  of  the  process  for  both  the  direct  method  and  the  first 


46 


guess  method  respectively.  The  result  from  the  first  guess  method  was  again  found  to  be 
even  more  accurate  than  the  direct  (single  assumption  set)  application.  A  linear 
regression  of  the  direct  method  returned  a  slope  of  1.3074,  an  offset  of  -153m  and  an 
overall  error  estimate  of  76m.  The  regression  of  the  first  guess  method  returned  a  line 
with  a  slope  of  1.2687,  an  offset  of -137m  and  an  overall  error  of  65m.  The  details  of 
this  application  and  how  it  compared  to  the  ground  truth  of  the  soundings  can  be  seen 
Table  2  in  Appendix  A. 


Actual  Boundary  Layer  Height  (m) 


Figure  4-14:  Linear  regression  of  actual  vs. 
calculated  boundary  layer  height  using  satellite 
derived  TB  data.  The  solid  line  is  the  regression  line. 
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Figure  4-15:  Linear  regression  of  actual  vs. 
calculated  boundary  layer  height  using  satellite  TB 
data.  The  solid  line  is  the  regression  line. 

*-  Heights  obtained  using  “first  guess”  process. 


1.  Favorable  Satellite  Case  Studies 

The  accuracy  of  using  satellite  Tb  as  cloud-top  temperature  in  the  technique 
described  above  will  be  illustrated  below  for  three  sounding  image  pairs. 

a.  GL16 

Figure  4-16  shows  the  closest  satellite  image  associated  with 
sounding  GL16  that  was  taken  at  2349  UTC  on  27  June  94.  The  visual  image  in 
the  area  of  this  sounding  shows  a  uniform  sheet  of  stratocumulus  cloud. 


48 


Figure  4-16:  AVHRR  Channel  2  image  for  GL16. 


Figure  4-17  shows  the  corresponding  TB  taken  from  the  IR  image  has  a  value  of 
1L4°C>  which  is  very  close  in  magnitude  to  the  sounding  cloud-top  temperature 
(1L0°C).  The  application  of  this  Tg  to  the  lapse  rate  calculation  of  BL  height 
yields  a  height  of  251m  (AH  =  -79m). 


49 


Figure  4-17:  AVHRR  channel  4  image  for  GL16. 

For  comparison  the  sounding  based  technique  had  a  height  differential  of 
AH  =  -30m  showing  that,  for  this  situation,  a  49m  change  in  accuracy  resulted  from 


a  Tct  error  of  +0.4°C. 


b.  GL67 

Figure  4-18  shows  the  image  that  most  closely  corresponds  to  sounding 
GL67  was  taken  at  1148  UTC  on  22Jun94.  The  corresponding  T3  taken  from  the  IR 
image  has  a  value  of  10.0°C,  which  is  very  close  in  magnitude  to  the  sounding  cloud-top 


temperature  (9.7°C). 
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Figure  4-18:  AVHRR  Channel  4  zoomed  image  from  for  GL67 


The  sounding  in  this  case  appears  to  go  through  a  region  of  thin  cloud  the  details  of  which 
can  be  clearly  seen  in  Figure  4-18.  The  corresponding  value  of  IB  measured  in  the  IR 
image  was  obtained  from  the  nearby  thicker  cloud  which  were  assumed  to  be 
representative  of  the  BL  in  the  vicinity  of  sounding  GL67.  The  application  of  this  Tb  to 
the  lapse  rate  calculation  of  BL  height  yields  a  height  of  477m  (AH  =  -29m).  For 
comparison  the  sounding  based  technique  had  a  height  differential  of  AH  =  -f3m  showing 
that,  for  sounding  GL67,  a  32m  change  in  accuracy  resulted  from  a  Tct  error  of  -K}.3°C. 
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GL68 


c. 

The  image  that  most  closely  corresponds  to  sounding  GL68  was  taken  at 
1514  UTC  on  22Jun94.  The  visual  image  (Figure  4-19)  in  the  area  of  this  sounding 
shows  a  region  of  even,  well  distributed  stratocumulus  clouds. 


Figure  4-19:  AVHRR  channel  2  image  for  GL68. 


The  corresponding  Tb  taken  from  the  IR  image  (Figure  4-20)  has  a  value  of  8.4°C.  which 
is  very  close  in  magnitude  to  the  sounding  cloud-top  temperature  (8.2°C).  The 
application  of  this  Tb  to  the  lapse  rate  calculation  of  BL  height  yields  a  height  of  647m 
(AH  =  -24m).  For  comparison  the  sounding  based  technique  had  a  height  differential  of 
AH  =  -6m  showing  that,  for  this  situation,  a  18m  change  in  accuracy  resulted  from  a  Tci 
error  of  +0.2°C. 
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Figure  4-20:  AVHRR  channel  4  image  for  GL68. 

Unfavorable  Satellite  Case  Studies 


There  are  several  situations  where  the  accurate  retrieval  of  Ter  from  satellite 
imagery  is  not  possible.  Errors  in  Tct  can  occur  when  there  are  higher  clouds  (e.g., 
cirrus)  above  the  STBL,  or  when  there  is  a  shallow'  BL  where  the  clouds  are  too  thin  and 
permit  radiance  from  the  surface  to  contaminate  the  TB  readings.  The  latter  case  will  be 
illustrated  below  as  will  an  instance  of  “limb  darkening”  where  TB  are  obtained  from  the 
far  edges  of  the  satellite  path. 
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a.  GL17 

The  closest  satellite  image  associated  with  sounding  GL17  was  taken  at 
1551  UTC  on  11  June  94.  This  image  can  be  seen  on  Figures  4-21  and  4-22.  which  show 
the  visual  (channel  2)  and  the  infrared  (channel  4)  respectively. 


Figure  4-21:  AVHRR  channel  2  image  for  GL1 7. 

The  brightness  temperature  derived  from  channel  4  for  this  image  was  10.8°C,  which  is 
slightly  warmer  than  the  cloud-top  temperature  measured  on  the  sounding  (10.1°C).  This 
increase  in  Tb  is  possible  due  to  the  influence  of  continental  aerosols  that  reduce  droplet 
size  and  absorptance  in  the  stratocumulus  clouds  leading  to  radiance  from  lower  in  the 
cloud  contributing  to  TB  and  producing  higher  values  of  Tg  (Brenner,  J.R.,  1994). 
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The  sounding  technique  results  that  were  described  previously  in  this 
chapter  returned  a  BL  height  of  286m  or  an  underestimate  of  AH  =  -64m.  The 
application  of  the  warmer  Tb  from  the  satellite  image  produced  an  enhancement  of  this 
underestimate  and  returned  a  BL  height  of  192m  or  a  AH  =  -157m.  This  indicates  that 
there  was  a  93m  change  in  accuracy  due  to  a  -K).7°C  change  in  Tct- 


Figure  4-22:  AVHRR  channel  4  image  for  GL17. 

b.  GL82 

Sounding  GL82  was  launched  at  2349  UTC  on  27  June  94  at  35.17°N- 
124.82°W.  The  surface  pressure  for  this  sounding  is  1021.8mb,  relative  humidity  is 
83.8%,  and  winds  are  from  the  NNW  at  9.2m/s.  The  surface  air/sea  temperature 
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difference  is  TA-TS  =  0.0°C.  and  the  surface  air  temperature  is  16.1°C.  The  depth  of  the 
boundary  layer  in  this  case  is  640m  with  a  cloud  fraction  of  approximately  34.3%. 

The  closest  satellite  image  associated  with  sounding  GL82  was  taken  at 
2334  UTC  on  27  June  94.  Figure  4-23  shows  the  DR.  satellite  image  for  GL82.  The  ima^e 
shows  that  the  corresponding  sounding  position  was  on  the  far  left-hand  side  of  the 
swath.  This  situation  falls  into  the  “limb  darkening”  category  where  radiance  from  the 
position  of  the  sounding  is  reduced  in  intensity  due  to  the  extreme  angle  at  which  it 
travels  to  reach  the  satellite  sensor. 


Figure  4-23:  AVHRR  channel  4  image  for  GL82. 

The  channel  4  brightness  temperature  measured  at  the  location  of  the 
sounding  on  Figure  4-23  has  a  value  of  9.1°C,  which  is  significantly  colder  than  the 
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corresponding  cloud-top  temperature  taken  from  sounding  GL82  (TCt  =  10.2°C).  This 
temperature  differential  results  in  a  calculated  BL  depth  of  808m  which  corresponds  to  a 
height  differential  of  AH  =  +168m. 

3.  Large  Scale  Satellite  Application 

The  next  step  in  the  logical  progression  of  this  technique  is  to  apply  it  to  the  entire 
image  vice  a  single  point.  This  was  done  using  a  single  Ts  measurement  from  sounding 
GL62  to  represent  the  entire  Ts  field  along  with  the  corresponding  AVHRR  image  shown 
in  Figures  2-24  and  4-25  that  were  taken  at  1535  UTC  on  21  June  94.  The  result  of  this 
process  can  be  seen  in  Figure  4-26  which  represents  a  topographic  map  of  BL  depth  off 
the  coast  of  California.  The  image  shows  the  depth  of  the  boundary  layer  increasing  away 
from  the  coast,  which  agrees  with  existing  theories  on  boundary  layer  depth  near 
coastlines  (Stull,  R.B.,  1988).  Depth  (in  meters)  can  be  determined  on  the  image  by  use 
of  the  color  scale  located  on  the  upper  left-hand  side  of  the  image.  The  image  contained 
higher  clouds  (seen  in  the  contrast  between  Figures  24  and  25)  and  cloud  free  areas  which 
were  filtered  out  and  appear  black  on  the  image, 
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Figure  4-24:  AVHRR  channel  2  image  for  sounding  GL62. 


Figure  4-25:  AVHRR  channel  4  Enhanced  IR  image  of  GL62. 
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Figure  4-26:  Topographic  map  of  BL  depth  based  on  sounding  GL62’s  SST, 
and  AVHRR  channel  4  imagery.  The  red  cross  indicates  the  location  of 
sounding  GL62. 
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V. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  purpose  of  this  thesis  was  to  develop  and  test  a  technique  for  determining  the 
depth  of  the  STBL  using  boundary  layer  thermodynamics.  The  technique  determines  BL 
depth  by  utilizing  satellite  derived  cloud-top  brightness  temperature  and  sea  surface 
temperature  observations.  The  method  was  initially  refined  with  sounding  data  from  the 
MAST  experiment  which  were  used  to  verify  the  technique.  The  technique  was  then 
applied  to  satellite  brightness  temperature  data  in  place  of  the  cloud-top  temperature  from 
the  soundings. 

Analysis  of  the  MAST  dataset  revealed  that  an  assumption  of  41%  cloud  (59% 
cloud-free)  with  a  moist  lapse  rate  of  Tm  =  -7.0°C/km  produced  the  best  overall  fit  to  the 
data.  However,  for  shallow  boundary  layers  the  most  accurate  assumption  set  was  75% 
cloud  with  a  lapse  rate  of  Tm  =  -6.5°C/km.  The  shallow  cases  were  found  to  apply  to 
boundary  layers  with  depths  less  than  400m.  The  value  of  400m  was  chosen  because  it 
was  found  to  be  the  center  of  a  gradual  shift  in  the  accuracy  of  the  two  assumption  sets. 

The  sensitivity  of  the  technique  was  analyzed  by  varying  the  values  of  both 
surface,  and  cloud-top  temperatures.  Analysis  revealed  that  a  14°C  change  in  either 
temperature  resulted  in  an  error  in  calculated  BL  depth  of  60-70m.  The  magnitude  of  BL 
height  error  was  found  to  vary  linearly  due  to  the  linearity  of  the  lapse  rates  used  in  the 
calculation  of  BL  depth.  The  sensitivity  of  this  technique  to  temperature  inputs  implies  a 
strong  need  for  accurate  Ts  and  TCT  fields  whenever  this  technique  is  used. 
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This  technique  is  applicable  to  regions  where  the  atmosphere  above  the  BL  is 
essentially  cloud  free  as  cirrus  clouds  significantly  impair  the  accurate  determination  of  TB 
from  the  STBL.  Furthermore,  the  technique  presented  here  was  designed  for  coupled  BLs 
in  a  uniform  cloud  field.  This  technique  will  not  yield  accurate  results  if  applied  to 
decoupled  BLs. 

B.  RECOMMENDATIONS 

This  study  quantitatively  demonstrates  the  sensitivity  of  BL  depth  determination  to 
the  accuracy  of  surface  and  cloud-top  temperature  data.  Models  that  generate  surface 
temperature  data  are  generally  of  much  lower  resolution  than  the  1km  accuracy  of 
AVHRR  Tb  data.  Future  work  in  this  area  should  ensure  that  the  resolution  of  both  Ts  and 
Tct  fields  are  compatible.  The  application  of  Ts  data  to  the  procedure  presented  in  this 
study  will  introduce  errors  due  to  differences  between  Ts  and  surface  air  temperature.  The 
errors  introduced  into  the  technique  due  to  this  difference  could  be  minimized  if  surface 
air  temperature  was  obtained  from  a  coupled  model  that  could  account  for  air/sea 
temperature  differences. 

Another  area  of  refinement  is  in  the  sensitivity  of  remotely  sensed  cloud-top 
temperature  to  cloud  thickness.  Measurements  of  TB  taken  from  thin  clouds  can  be 
contaminated  with  surface  radiance,  which  makes  the  value  of  TB  too  high.  Multi-channel 
remote  sensing  techniques  can  be  employed  to  estimate  cloud  optical  thickness  and  correct 
the  measured  TB  to  account  for  contamination  from  radiance  below  the  cloud-top. 
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The  application  of  this  technique  to  EM  propagation  and  ducting  conditions  can  be 
refined  with  the  application  of  multi-spectral  remote  sensing  techniques.  Multi-spectral 
techniques  could  be  employed  to  enhance  the  determination  of  the  moisture  content  above 
the  boundary  layer  which  will  help  with  the  determination  of  refractivity  conditions  (i.e., 
ducting). 

Additional  study  should  be  done  with  this  technique  in  order  to  better  determine 
the  nature  of  shallow  BLs  (<400m)  and  BLs  deeper  than  700m.  Boundary  layers 
shallower  and  deeper  than  400m  and  700m  respectively  were  not  well  represented  in  this 
study.  The  technique  presented  in  this  study  was  not  sufficiently  sophisticated  to  deal 
with  complexity  inherent  in  decoupled  boundaiy  layers.  Therefore,  another  area  for 
further  research  is  the  development  of  a  technique  to  account  for  decoupling  processes 
within  the  STBL. 

The  ability  to  estimate  boundary  layer  depth  from  satellite  data  is  a  valuable 
additional  tool  for  mapping  the  coastal  battlespace  or  operations  areas.  In  addition, 
boundary  layer  depth  analyses  could  provide  valuable  additional  data  for  local  mesoscale 
models  that  have  been  deployed  recently.  This  study  points  to  the  way  for  an  important 
new  qualitative  use  of  satellite  data  to  support  Navy  operations. 
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APPENDIX  A.  TABLES 


This  appendix  serves  as  a  convenient  location  in  order  to  consolidate  the  tables 
discussed  in  the  text  of  this  thesis. 
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Table  1:  Cloud  amount  and  lapse  rate  assumptions  and  how  they  affect  the 
accuracy  when  compared  to  the  actual  height  of  the  STBL. 


Cloud  Amount 
(fraction  of  BL) 

Lapse  Rate 
(°C/Km) 

Linear 

Regression 

Slope 

Y-intercept 
offset  (m) 

Overall  Error 
Estimate  (m) 

75%  Cloud 
<400m) 

41%  Cloud 
>400m) 

-6.5  <400m 
-7.0  >400m 

1.0283 

-41 

56 

75%  Cloud 
<400irT 

41%  Cloud 
>400m* 

-6.5  <400m  * 
-7.0  >400m* 

1.0080 

-24 

50 

75%  Cloud 

-7.0 

1.2174 

-93 

60 

75%  Cloud 

-6.5 

1.2949 

-103 

64 

60%  Cloud 

-7.0 

1.1488 

-86 

58 

60%  Cloud 

-6.5 

1.2117 

-98 

61 

41%  Cloud 

-7.0 

1.0861 

-82 

53 

41%  Cloud 

-6.5 

1.1236 

-85  | 

54 

*-  Denotes  ’’first  guess”  process  to  determine  the  400m  dividing  line  for  the  two 
sets  of  assumptions. 
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Table  2:  A  comparison  of  how  the  accuracy  of  the  technique  was  affected  by 
using  Tb  for  the  TCT  vice  sounding  based  Tct- 


Cloud 
Amount 
(fraction  of 
boundary 
layer) 

Lapse  Rate 
CC/Km) 

Linear 

Regression 

Slope 

Y- 

intercept 
offset  (m) 

Overall 

Error 

Estimate 

(m) 

Sounding 

Based 

Results** 

All 

Soundings 

75%Cloud 
(<400m) 
41%  Cloud 
(>400m) 

-6.5  (<400m) 
-7.0  (>400m) 

1.0080 

-24 

50 

Sounding 

Based 

Results** 

Satellite 

Equivalent 

Sample 

Same  as 
above 

Same  as 
above 

0.9975 

-23 

49 

Satellite  TB 
Results* 

Same  as 
above 

Same  as 
above 

1.3074 

-153 

76 

Primary 
Satellite  TB 
Results** 

Same  as 
above 

Same  as 
above 

1.2687 

-137 

65 

The  “satellite  equivalent  sample”  indicates  that  the  same  soundings  used  for  the 
satellite-based  results  were  extracted  from  the  full  sounding  set  and  run 
separately. 

*  -  Indicates  the  assumptions  were  divided  according  to  the  previously  known 
(actual)  height  of  the  BL. 

**  -  Denotes  ’’first  guess”  process  to  determine  the  400m  dividing  line  for  the 
two  sets  of  assumptions. 
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Table  3:  MAST  experiment  soundings,  with  their  corresponding  position, 
time  and  analysis. 


Sounding 

Number 

Latitude, 

Longitude 

Launch  Time 
(date  -  UTC) 

Analysis 

GL01 

34.15°N 

119.21°W 

03Jun  -  0029 

Dismissed  -  test 
sounding 

GL02 

34.16°N 

119.27°W 

05Jun  -  2357 

Dismissed  -  no  cloud 

GL03 

35.25°N 

121.16°W 

06Jun  - 1447 

Dismissed  -  no  cloud 

GL04 

35.98°N 

121.65°W 

-2039 

Dismissed  -  no  cloud 

GL05 

36.63°N 

122.00°W 

08Jun  -  0535 

Dismissed  -  no  cloud 

GL06 

35.73°N 

122.22°W 

- 1149 

Dismissed  -  no  cloud 

GL07 

35.16°N 

122.44°W 

- 1737 

Dismissed  -  no  cloud 

GL08 

35.29°N 

122.46°W 

-2353 

Dismissed  -  no  cloud 

GL09 

36.61°N 

121.89°W 

lOJun-  1739 

Dismissed  -  no  cloud 

GL10 

36.45°N 

122.07°W 

- 1929 

Dismissed  -  no  cloud 

GL11 

36.27°N 

122.18°W 

-2058 

Dismissed  -  no  cloud 

GL12 

36.08°N 

122.18°W 

-2354 

Conducive 

GL13 

36.25°N 

122.48°W 

llJun-0252 

Dismissed  -  too  shallow 

GL14 

36.26°N 

122.99°W 

-0549 

Conducive 

GL15 

36.26°N 

123.00°W 

-0838 

Dismissed  -  bad 
sounding 

GL16 

36.15°N 

122.57°W 

- 1153 

Conducive 

GL17 

36.23°N 

122.48°W 

- 1448 

Conducive 

GL18 

36.52°N 

122.30°W 

- 1746 

Conducive 
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Table  3  continued. 


GL19 

36.35°N 
122.3 1°W 

-2041 

Conducive 

GL20 

36.66°N 
122.4 1°W 

-2343 

Conducive 

GL21 

36.86°N 

122.67°W 

12Jun  -  0244 

Dismissed  -  no  cloud 

GL22 

37.07°N 

122.97°W 

-0541 

Dismissed  -  no  cloud 

GL23 

37.3°N 

123.26°W 

- 1146 

Dismissed  -  bad 
sounding 

GL24 

37.23 °N 
123.22°W 

- 1448 

Dismissed  -  bad 
sounding 

GL25 

36.99°N 

122.89°W 

- 1732 

Conducive 

GL26 

36.90°N 

123.06°W 

-2020 

Conducive 

GL27 

36.96°N 

123.14°W 

-2132 

Conducive 

GL28 

37.05°N 

123.25°W 

-2347 

Conducive 

GL29 

37.19°N 

123.40°W 

13Jun  -  0242 

Conducive 

GL30 

36.95°N 

123.01°W 

-0549 

Conducive 

GL31 

36.65°N 

122.77°W 

-0841 

Conducive 

GL32 

36.3 1°N 
122.69°W 

- 1143 

Conducive 

36.01°N 

122.60°W 

- 1451 

Conducive 

GL34 

36.07°N 

122.56°W 

- 1741 

Dismissed  -  bad 
sounding 

GL35 

36.05°N 

122.75°W 

-2030 

Dismissed  -  no  cloud 

GL36 

36.33°N 

123.16°W 

-2338 

Dismissed  -  no  cloud 

GL37 

36.35°N 

122.65°W 

14Jun  -  0242 

Dismissed  -  no  cloud 

GL38 

36.53°N 

122.16°W 

-0537 

Dismissed  -  no  cloud 
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Table  3  continued. 


GL39 

36.14°N 

122.08°W 

15Jun-  1740 

Dismissed  -  no  cloud 

GUO 

35.46°N 

122.26°W 

-2324 

Dismissed  -  no  cloud 

GUI 

35.50°N 

122.42°W 

16Jun  -  0257 

Dismissed  -  no  cloud 

GU2 

35.69°N 

122.61°W 

-0542 

Dismissed  -  no  cloud 

GU3 

35.20°N 

123.25°W 

-  1145 

Conducive 

GU4 

35.32°N 

123.64°W 

- 1456 

Dismissed  -  no  cloud 

GU5 

35.50°N 

123.78°W 

-  1735 

Dismissed  -  no  cloud 

GU6 

35.6 1°N 
123.90°W 

-2030 

Dismissed  -  no  cloud 

GU7 

35.85°N 

124.00°W 

-2335 

Dismissed  -  no  cloud 

GU8 

36.09°N 

124.12°W 

3 

o 

1 

o 

Dismissed  -  no  cloud 

GU9 

36.30°N 

124.23°W 

-0545 

Dismissed  -  no  cloud 

GL50 

36.5 1°N 
122.33°W 

-  0841 

Dismissed  -  no  cloud 

GL51 

36.73°N 

124.41°W 

- 1156 

Conducive 

GL52 

36.85°N 

124.29°W 

-  1448 

Conducive 

GL53 

36.85°N 

123.76°W 

-  1804 

Dismissed  -  no  cloud 

GL54 

36.76°N 

122.73°W 

-2348 

Dismissed  -  no  cloud 

GL55 

36.72°N 

122.01°W 

20Jun  -  1458 

Conducive 

GL56 

36.77°N 

122.42°W 

-  1809 

Conducive 

GL57 

36.56°N 

122.59°W 

-2358 

Dismissed  -  no  cloud 
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Table  3  continued. 


GL58 

36.24°N 

122.54°W 

21Jun  -  0233 

Dismissed  -  no  cloud 

GL59 

35.88°N 

122.37°W 

-0551 

Dismissed  -  no  cloud 

GL60 

35.54°N 

122.13°W 

-0847 

Conducive 

GL61 

35.15°N 

121.90°W 

- 1150 

Conducive 

GL62 

35.1 1°N 
121.85°W 

- 1503 

Conducive 

GL63 

35.1 1°N 
121.88°W 

- 1751 

Conducive 

GL64 

34.83°N 

121.77°W 

-2053 

Conducive 

GL65 

35.02°N 

1 22.06°  W 

-2352 

Dismissed  -  no  cloud 

GL66 

35.22°N 

122.27°W 

22Jun  -  0244 

Dismissed  -  no  cloud 

GL67 

35.78°N 

123.18°W 

-  1148 

Conducive 

GL68 

35.70°N 

123.46°W 

-  1453 

Conducive 

GL69 

35.26°N 

123.58°W 

-  1752 

Conducive 

GL70 

34.98°N 

123.73°W 

-2046 

Dismissed  -  no  cloud 

GL71 

35.27°N 

123.65°W 

-2346 

Dismissed  -  no  cloud 

GL72 

35.58°N 

123.32°W 

23  Jun  -  0249 

Dismissed  -  no  cloud 

GL73 

35.88°N 

123.02°W 

-0545 

Dismissed  -  no  cloud 

GL74 

35.50°N 

1 22.06 °W 

-  2006 

Dismissed  -  no  cloud 

GL75 

36.65°N 

121.94°W 

26Jun  -  0546 

Dismissed  -  no  cloud 

GL76 

35.90°N 

123.01°W 

-  1453 

Dismissed  -  no  cloud 
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Table  3  continued. 


GL77 

35.10°N 

123.67°W 

-2151 

Dismissed  -  no  cloud 

GL78 

34.53°N 

124.02°W 

27Jun-  0241 

Conducive 

GL79 

34.40°N 

124.78°W 

-1452 

Conducive 

GL80 

34.54°N 

124.84°W 

-1801 

Conducive 

GL81 

34.78°N 

124.83°W 

-2053 

Conducive 

GL82 

35.17°N 

124.84°W 

-2349 

Conducive 

GL83 

35.44°N 

124.82°W 

28Jun  -  0238 

Conducive 

GL84 

35.71°N 

124.66°W 

-0545 

Conducive 

GL85 

35.33°N 

124.52°W 

-0855 

Conducive 

GL86 

35.03°N 

124.03°W 

-1150 

Conducive 

GL87 

35.1 1°N 
123.89°W 

-  1451 

Conducive 

GL88 

35.33°N 

124.07°W 

-  1800 

Conducive 

GL89 

35.55°N 

124.25°W 

-2057 

Conducive 

GL90 

35.74°N 

1 24.44°  W 

-2344 

Conducive 

GL91 

35.91°N 

124.58°W 

29Jun  -  0249 

Conducive 

GL92 

35.58°N 
124.3 7°W 

-0547 

Conducive 

GL93 

35.38°N 

124.00°W 

-0841 

Conducive 

GL94 

34.94°N 
122.7 1°W 

-  1645 

Conducive 
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Table  4:  MAST  experiment  soundings,  corresponding  satellite  images  and 
analysis. 


Sounding 

Number 

Launch  Time 
(date  -  UTC) 

Corresponding 
Satellite  Image 
(date;  sat;  UTC) 

Analysis 

GL12 

lOJun  -  2354 

llJun;  nil;  0119 

Conducive 

GL14 

1 1  Jun  -  0549 

n9;  0530 

Conducive 

GL16 

- 1153 

nil; 1342 

Conducive 

GL17 

- 1448 

nl2; 1551 

Conducive 

GL18 

- 1746 

n9; 1753 

Conducive 

GL19 

-2041 

none 

Dismissed  -  no  image 

GL20 

-2343 

nil;  2327 

Conducive 

GL25 

12Jun-  1732 

12 Jun;  n9;  1740 

Dismissed  -  cirrus 
contamination 

GL26 

-2020 

none 

Dismissed  -  no  image 

GL27 

-2132 

none 

Dismissed  -  no  image 

GL28 

-2347 

13Jun;  nil;  0057 

Dismissed  -  cirrus 
contamination 

GL29 

13  Jun  -  0242 

nl2;  0249 

Dismissed  -  cirrus 
contamination 

GL30 

-0549 

n9;  0503 

Dismissed  -  cirrus 
contamination 

GL31 

-0841 

none 

Dismissed—  no  image 

GL32 

-  1143 

nil; 1318 

Dismissed  -  cirrus 
contamination 

GL33 

-  1451 

nl2;  1508 

Dismissed  -  cirrus 
contamination 

GL43 

16Jun- 1145 

16Jun;  nil;  1241 

Dismissed  -  cirrus 
contamination 

GL51 

17Jun- 1156 

17Jun;  nil;  1229 

Dismissed  -  cirrus 
contamination 

GL52 

- 1448 

nil; 1410 

Dismissed  -  cirrus 
contamination 
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Table  4  continued. 


GL55 

20Jun  -  1458 

20Jun;  nlO;  1455 

Dismissed  -  cirrus 
contamination 

GL56 

- 1809 

n9; 1737 

Conducive 

GL60 

21Jun  -  0847 

none 

Dismissed  -  no  image 

GL61 

- 1150 

21Jun;  nil;  1320 

Conducive 

GL62 

-  1503 

N12;  1535 

Conducive 

GL63 

- 1751 

N9;  1724 

Conducive 

GL64 

-2053 

None 

Dismissed  -  no  image 

GL67 

22Jun  -  1 148 

22Jun;  nil;  1308 

Conducive 

GL68 

-  1453 

nl2; 1514 

Conducive 

GL69 

- 1752 

n9; 1711 

Conducive 

GL78 

27Jun-  0241 

27Jun;  nl2;  0246 

Dismissed  -  cirrus 
contamination 

GL79 

-  1452 

nl2; 1506 

Dismissed  -  limb 
darkened 

GL80 

-  1801 

n9; 1747 

Conducive 

GL81 

-.2053 

none 

Dismissed  -  no  image 

GL82 

-2349 

nil;  2332 

Dismissed  -  limb 
darkened 

|| 

28Jun  -  0238 

28Jun;  nl2;  0246 

Conducive 

-0545 

n9;  0511 

Conducive 

GL85 

-0855 

none 

Dismissed  -  no  image 

GL86 

- 1150 

none 

Dismissed  -  no  image 

GL87 

nlO;  1503 

Conducive 

GL88 

n9;  1734 

Conducive 

GL89 

-2057 

none 

Dismissed  -  no  image 

GL90 

-2344 

nil;  2320 

Dismissed  -  no  data 

29Jun  -  0249 

nlO;  0221 

Conducive 

-0547 

29Jun;  n9;  0458 

Conducive 

GL93 

-0841 

none 

Dismissed  -  no  image 

GL94 

-  1645 

n9; 1721 

Conducive 
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APPENDIX  B.  FORMULAE  AND  PROCEDURES 


The  following  description  of  the  formulae  and  procedures  used  in  processing  the 
sounding  data  from  the  MAST  experiment  was  provided  by  William  Syrett  of 
Pennsylvania  State  University. 

The  measured  parameters  were  pressure,  temperature  and  relative  humidity,  with 
winds  calculated  using  the  Omega  network.  The  "Omega"  winds  were  averaged  over  a 
four-minute  period  by  the  sounding  system  software.  Raw  temperature  and  humidity  data 
were  output  at  1.5  second  intervals,  with  winds  output  every  10  seconds.  The  raw 
temperatures  and  humidities  were  cleaned  up  and  then  averaged  to  5-second  intervals, 
while  winds  were  simply  interpolated  to  5  seconds. 

Derived  quantities  include  height,  dewpoint  and  mixing  ratio.  The  formulas  and 
procedures  used  for  mixing  ratio,  dewpoint  and  height  are  (in  FORTRAN  format): 

Mixing  Ratio  (w): 

w  =  622.0*(e/(p-e))  p  =  pressure  (mb) 

e  =  vapor  pressure  (mb) 

where:  e  =  RH*es/100.0  RH=  relative  humidity 

es=  sat.  vapor  pres,  (mb) 

where:  es=  6.1 12*EXP((17.67*T)/(T+243.5)) 

T  =  temperature  (deg  C) 

Height  (z): 

z(i)  =  z(i- 1  )+(R*T vA/ g)  *LN (p(i- 1  )/p(i))  R  =  gas  constant  ( J/kg  K) 

g  =  acc.  due  to  gravity 

where:  TvA  is  the  layer-averaged  virtual  temp; 
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Tv  =  T*(l. 0+0.6  l*w) ;  w  in  g/g,  Tv  in  degrees  Kelvin  for  height 
computation. 

Dewpoint: 

Td  =  (243 . 5  *LN  (e/6. 1 1 2))/(  1 7 . 67-LN(e/6 . 1 1 2)) 

Please  note  that  the  relative  humidities  were  adjusted  upward,  based  on 
both  observed  cloudiness  in  relation  to  reported  humidities  and  also  on  a  conversation  with 
a  Vaisala  employee  familiar  with  the  humidity  sensor  on  the  RS-80  sondes  used  at  the 
sites.  The  adjustment  procedure  is  similar  to  that  used  for  ASTEX  soundings,  except  the 
maximum  upward  adjustment  has  been  increased  to  7%,  up  from  5%.  The  adjustment 
procedure  is  described  next. 

***  Upward  Adjustment  of  Relative  Humidity  Values  *** 

An  initial  data  quality  check  was  run.  A  maximum  value  was  determined  from 
the  accepted  relative  humidities.  If  the  maximum  unadjusted  humidity  was  greater  than  or 
equal  to  65%  a  correction  factor  was  added  to  the  maximum  humidity.  This  factor  ranges 
from  1%  if  the  maximum  unadjusted  humidity  was  65%,  to  7%  if  the  highest  humidity 
was  93%.  Maxima  at  or  above  93%  were  adjusted  to  100%. 

This  correction  factor  (C.F.)  was  then  linearly  weighted  according  to  the  humidity. 
The  equation  for  humidity  adjustment  is: 

RH(new)  =  RH(old)  +  ((RH(old)  -  20.0)/(RH(max)  -  20.0))  *  C.F. 

There  was  no  adjustment  if  the  humidity  was  less  than  20%. 

Obvious,  singular  bad  data  (spikes  in  humidity  and  temperature)  were  removed. 
Unstable  lapse  rates  found  in  the  lowest  layers  were  the  result  of  insufficient  exposure  of 
the  sonde  to  the  ambient  environment  before  launch  (primarily  early  in  the  campaign)  and 
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were  thus  removed.  If  the  surface  temperature  had  to  be  lowered  more  than  0.5  degrees 
Celsius  a  comment  was  generated  in  the  summary  table.  Only  the  lowest  levels  were 
adjusted  as  by  20  seconds  into  launch  the  temperature  sensor  generally  had  "caught  up" 
with  the  actual  (generally  dry-adiabatic)  profile. 

The  general  data  quality  was  quite  good.  Temperature  data  were  consistently  of 
high  quality  with  the  exception  of  the  false  super-adiabatic  layers.  Humidity  data  were 
generally  of  good  quality,  although  values  were  a  few  percent  low  when  the  air  was  moist. 
The  humidity  profiles  for  several  of  the  first  approximately  35  soundings  appeared 
"unrealistic"  in  real  time.  The  problem,  if  it  is  one,  manifested  itself  as  regions  of  higher 
relative  humidity  just  above  the  inversion.  The  resulting  mixing  ratio  profiles  appeared 
unrealistic.  Specifically,  soundings  10, 1 1, 23, 32, 33  and  34  are  suspect-  use  with 
CAUTION! 

Note:  Soundings  10,  11,  23,  32,  33,  and  34  were  analyzed  more  closely  and 
soundings  32  and  33  were  deemed  useful  for  the  purposes  of  this  study. 

Low-level  winds  were  unfortunately  difficult  to  obtain.  It  was  often  difficult  to 
get  a  "lock",  necessary  for  low-level  winds,  while  the  sonde  was  on  board.  Winds 
generally  began  being  received  when  the  sonde  was  above  1  km  in  altitude. 
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